function [d, t] = getcmg(dat, th, cmpNo)
mp = arrayfun(@(x) (x.SourceX+x.GroupX)/2, th);
d = dat(:, mp==cmpNo);
t = th(mp==cmpNo);
end